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Abstract (Basic; : EP 1074858 Al 

NOVELTY - The pseudo -period corresponding to the inverse of the 
dominant frequency of the frequency spectrum is determined and binary 
coded relative to a first threshold value (THl) . 

The coded signals are summed through a sliding window to give an 
integral signal (S) and a second threshold value (TH2) . The significant 
instant is located when the corresponding integral signal reaches the 
second threshold value 

DETAILED DESCRIPTION - The location of "the significant instant of " 
each signal involves determining a first reference point where the 
amplitude of the corresponding integral signal is equal to the second 
threshold value. A first intermediate instant (LP0) is detected between 
- this first reference point and the nearest milestone on the integral 
signal. A second intermediate instant is detected between the first 
intermediate instant and the opposite milestone of the integral signal, 
where a first maximum (SM) of this integral signal exceeds a third 
threshold value (TH3) depending on the acquired signal. The significant 
instant is located at a quarter pseudo-period from the second 
intermediate instant. The significant instant is the first arrival of 
each acquired signal or the end of each acquired signal. The signals 
between these instants are isolated on the recording of the data. Each 
part of the data between two such instants is isolated so all the 
separate arrivals are identified for each detector. The polarisation of 
the signals is also detected so compression waves, sent waves and 
multiplets can be identified. 

USE - Analysis of signals received by one detector from a number of 
sources or by several detectors from one source, for example in seismic 
surveys 

ADVANTAGE - Accurate- detection of significant instants allows their 
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source to be located. Current methods involve the operator identifying- 
the instants manually in a long and fastidious process 
pp; 13 DwgNo 0/5 
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Description 

[0001] La presente invention concerne une methode 
d'analyse de signaux regus sur des capteurs d'ondes 
elastiques permettant par exemple de pointer automa- 5 
tiquement avec une grande precision des instants signi- 
ficatifs tels que les instants de premiere arrivee sur les 
differents capteurs et/ou les instants de fin de ces me- 
rries signaux. 

[0002] La methode peut trouver des applications dans 
de nombreux domaines ou Ton doit pointer avec une 
grande precision les instants d'arrivee de signaux vibra- 
toires. C'est le cas notamment dans le cadre de syste- 
mes de surveillance sismique active ou passive de gi- 
sements souterrains en cours d'exploitation utilises 
d'une facon generale soit pour en extraire des fluides, 
soit pour y injecter des fluides. 

ETAT DE LA TECHNIQUE 

[0003] Par les brevets FR 2 593 292 (US 4,775,009), 
FR 2 681 373 (US 5,303,773), FR 2 685 139 (US 
5,370,546) et FR 2 703 470 notamment, on connait dif- 
ferentes techniques pour surveiller revolution au cours 
du temps de reservoirs souterrains comportant I'utilisa- 
tion de capteurs sismiques ou autres installes a demeu- 
re dans un ou plusieurs puits, sans perturber les opera- 
tions diverses (production, ou injection de fluides, inter- 
ventions diverses menees par I'intermediaire de ces 
puits). L'installation de ces capteurs a demeure dans 
des puits, permet d'exercer une surveillance sismique 
d'un reservoir dans le but de detecter differents pheno- 
menes lies a son exploitation. Les capteurs sismiques 
permanents sont installes par exemple a I'exterieur d'un 
casing ou cuvelage que Ton descend dans le puits. Ms 
sont noyes dans le ciment qui est normalement injecte 
dans I'espace annulaire entre lui et le puits et qui assure 
leur bon couplage avec les formations environnantes. 
Les capteurs sismiques peuvent encore etre fixes a I'ex- 
terieur d'une colonne d'exploitation qui est installee 
dans un putts lors des operations de completion, lis sont 
associes a des moyens mobiles adaptes a les plaquer 
contre le casing du puits, et a des moyens de decoupla- 
ge pour fiitrer les vibrations se propageant le long de la 
colonne d'exploitation. 

[0004] Par les brevets FR 2 703 457, FR 2 703 470 
ou FR 2 728 973 (US 5,724,311 ), on connait aussi des 
methodes de surveillance active repetitive a long terme 
d'un reservoir par application d'ondes elastiques a une 
formation et acquisition des signaux-reponses renvoy6s 
par la formation, ceci au moyen d'une installation per- 
manent^ de moyens d'emission et de reception dans 
des puits ou au voisinage de la surface du sol. On realise 
des traitements differentiels sur des acquisitions op6- 
rees dans des conditions identiques. 
[0005] Par le brevet FR 2 772 137, on connaTt aussi 
une methode permettant de simplifier les problemes de 
discrimination et d'identification, d'evenements sismi- 



ques ou microsismiques relatifs a I'exploitation d'une zo- 
ne souterraine. On utilise un ou plusieurs capteurs de 
reference presentant un couplage acoustique direct 
avec des elements de I'equipement technique d'exploi- 
tation de la zone qui peuvent detecter aussi des ondes 
liees directement a I'exploitation. L'analyse comparee 
des signaux issus des differents recepteurs et de cha- 
que capteur de reference, permet de trier les enregis- 
trements en differentes families selon que les evene- 
ments dans la zone souterraine sont independants des 
evenements detectes par chaque capteur de reference 
ou en dependent directement ou indirectement. On peut 
etablir un catalogue d'evenements-types ayant un lien 
causal avec des interventions d'exploitation sur la zone 
souterraine, par correlation entre les enregistrements 
des signaux regus respectivement par les dits recep- 
teurs et le (ou chaque) capteur de reference, et ainsi 
identifier la nature des evenements sismiques ou micro- 
sismiques. 

[0006] La determination precise des instants d'arrivee 
des signaux regus par plusieurs recepteurs couples 
avec la formation souterraine rend possible la localisa- 
tion des lieux d'emission (points-sources ou foyers) ou 
se produisent des evenements qui sont lies directement 
ou indirectement a une activite ayant pour effet de mo- 
difier le champ de contraintes : operations d'injeqtion ou 
de soutirage de fluides par un ou plusieurs puits par 
exemple. 

[0007] Ceci permet a I'exploitant d'un site d'interpreter 
tres rapidement les donnees de position obtenues en 
relation avec les donnees classiques d'exploitation 
(pression, debit, temperature etc.) de sorte qu'il peut 
prendre en compte la reponse mecanique du site dans 
le protocole d'exploitation afin de preserver la producti- 
vity du ou des puits voire du site. 
[0008] Cependant, la methode la plus courante pour 
detecter avec precision les instants significatifs des si- 
gnaux regus lors d'operations de surveillance sismique 
d'une zone souterraine consiste essentiellement pour 
I'operateur a visionner tous les enregistrements obte- 
nus et a pointer lui-meme les instants (de premiere ar- 
rivee) ou ils se produisent sur tous les capteurs. C'est 
la une operation qui s'avere longue et fastidieuse quand 
on doit depouiller une grande masse d'enregistrements 
systematiques effectues sur de longues dur6es. 

LA METHODE SELON ^INVENTION 

[0009] La methode selon I'invention d'analyse de si- 
gnaux acquis provenant d'au moins un capteur d'ondes, 
permet de detecter automatiquement avec une grande 
precision, au moins un instant significatif de chaque si- 
gnal tel que I'instant de sa premiere arrivee (ti) ou I'ins- 
tant de sa fin (tf), et done de simplifier les operations de 
depouillement et de tri prealables a Interpretation. 
[0010] Elle comporte : 

la determination de la pseudo-periode correspon- 
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dant a I'inverse de la frequence dominante du spec- 
tre de frequence d'au moins une partie de chaque 
signal acquis, 

le codage de chaque signal acquis pour former un 
signal code en binaire par reference a une premiere 
valeur-seuil dependant du signal acquis (codage a 
la valeur 1 ou 0 par exemple, selon que sa valeur 
absolue est au-dessus ou au-dessous de cette va- 
leur-seuil), 

la formation pour chacun d'eux d'un signal integral 
par sommation du signal code en binaire sur une 
fenetre glissante de maniere a determiner une 
deuxieme valeur-seuil (dependant de preference 
de I'amplitude maximale de ce signal integral), et 

ie positionnement de I'instant significatif sur le si- 
gnal acquis par reference a un instant ou le signal 
integral correspondant atteint cette deuxieme va- 
leur-seuil. 

[001 1] Pour positionner I'instant significatif de chaque 
signal acquis on determine par exemple un premier 
point de reference ou I'amplitude du signal integral cor- 
respondant est egal a cette deuxieme valeur-seuil, on 
detecte un premier instant intermediate entre ce pre- 
mier point de reference et la borne la plus proche du dit 
signal integral (le debut, si Ton cherche a detecter I'ins- 
tant de premiere arrivee), on detecte aussi un deuxieme 
instant intermediaire entre ce premier instant interme- 
diaire et la borne opposee du signal integral (la fin, dans 
la meme hypothese) ou ce signal integral depasse une 
troisieme valeur-seuil dependant du signal acquis. On 
repere alors I'instant significatif recherche sur le signal 
acquis a un quart de pseudo-periode du deuxieme ins- 
tant intermediaire (en amont ou en aval selon qu'il s'agit 
de I'instant de la premiere arrivee ou de I'instant de fin 
de chaque signal acquis). 

[0012] On peut a la fois determiner I'instant de la pre- 
miere arrivee de chaque signal acquis et I'instant de fin 
de chaque signal acquis et isoler sur les enregistre- 
ments des signaux regus par chaque capteur d'ondes, 
des parties d'enregistrement entre ces deux instants. 
[0013] On peut aussi, sur un meme signal acquis, iso- 
ler toutes les parties d'enregistrement delimitees par dif- 
ferents instants significatifs de maniere a detecter tou- 
tes les arrivees distinctes sur un meme capteur d'ondes. 
[001 4] La methode peut comporter aussi une analyse 
des signaux regus respectivement par une plurality de 
capteurs d'ondes relatifs a un meme ev6nement gene- 
rates de signaux elastiques, une determination de la 
polarisation de chaque signal acquis immediatement 
apres son instant de premiere arrivee et une detection 
des arrivees successives d'ondes de compression et 
d'ondes de cisaillement par comparaison des orienta- 
tions des vecteurs propres lies aux polarisations respec- 
tives de ces signaux. 



[001 5] La methode peut aussi comporter une analyse 
des signaux acquis provenant respectivement d'une 
plurality de capteurs d'ondes relatifs a un meme evene- 
ment generateur de ces signaux, une determination de 
5 la polarisation de chaque signal acquis immediatement 
apres son instant de premiere arrivee et une detection 
de multiplets par comparaison des orientations des vec- 
teurs propres lies aux polarisations respectives de ces 
signaux. 

10 [0016] La methode selon I'invention permet la detec- 
tion automatique de I'arrivee d'une onde (notamment 
sismique) avec une tres bonne precision. Elle est done 
auto-adaptative car elle atteint son objectif : evaluation 
de la frequence dominante de chaque signal acquis et 

15 utilisation de la pseudo-periode correspondant a cette 
frequence dominante pour investiguer un signal integral 
(code en binaire et somme) etabli a partir du signal initial 
brut, en tenant compte pour cela d'un jeu de quelques 
parametres ajustables en fonction des conditions loca- 

20 les d'acquisition des signaux, mais dependant du signal 
lui-meme. 

[0017] Elle facilite done la localisation automatique 
des evenements microsismiques presentant une suc- 
cession d'ondes de compression et de cisaillement et, 

25 si ces evenements sont detectes par un nombre suffi- 
sant de capteurs espaces elle permet une triangulation 
automatique a partir des temps d'arrivee ainsi determi- 
nes, le modele de vitesse etant connu par ailleurs. 
[0018] D'autres caracteristiques et avantages de la 

30 methode selon I'invention, apparaitront a la lecture de 
la description ci-apres d'un exemple non limitatif de rea- 
lisation, en se referant aux dessins annexes ou : 

les Fig. 1a a 1c montrent un exemple de signaux 
35 captes respectivement sur les trois composantes 
d'un geophone a trois voies ( dit triphone) ; 

les Fig.2a a 2c montrent respectivement les si- 
gnaux de la Fig.1 a plus grande echelle ; 

40 

la Fig.3a montre le codage binaire, applique a un 
signal, a 1 ou 0 selon qu'il depasse ou non un ni- 
veau-seuil ; 

45 - la Fig. 3b montre le signal integral forme par la som- 
mation du signal code de la Fig.3a, sur une fenetre 
glissante ; 

les Fig. 4a a 4c montrent respectivement les si- 
50 gnaux integraux formes a partir des signaux res- 
pectivement des Fig. 1a a 1c ou 2a a 2c ; et 
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DESCRIPTION DETAILLEE 

I) Etapes preiiminaires 
[0019] 

1-1) On acquiert numeriquement chaque signal 
avec une frequence d'echantillonnage suffisante 
f ech (2400 Hz par exemple dans I'exemple pr6sent6) 
depuis un instant initial pr6c6dant suffisamment les 
possibles premieres arrivees de faible amplitude de 
facon a ne pas les manquer (zone de pr6-deden- 
chement). 

I-2) De preference, on filtre le signal acquis par cha- 
que capteur sans le dephaser pour ne conserver 
que le spectre de frequence dont on aura au pr£a- 
lable fix£ les limites superieure et inferieure "fmin" 
et "fmax 0 (parametres modifiables fixes par exem- 
ple a 5 Hz et 300 Hz). 

I- 3) De preference, on filtre egalement chaque si- 
gnal acquis au moyen d'un filtre median param6tr6 
MF afin d'eiiminer les pics parasites sortant du ni- 
veau de bruit et de dur£e voisine de la periode 
d'echantillonnage du signal. 

Les deux operations de filtrage precedentes 
sont recommandees mais facultatives. 

1 -4) On calcule la frequence dominante F dom dans 
le signal par une technique connue d'analyse spec- 
trale (FFT par exemple). S'il y a des ecarts entre les 
valeurs de F dom obtenues sur les differentes voies 
d'un capteur a trois composantes, on determine une 
valeur moyenne. On en deduit la pseudo-periode 
correspondante : T PS =1/ F dom , ce qui permet de 
determiner le nombre de periodes (points) d'6chan- 
tillonnage np par pseudo-p6riode T PS de signal : 
"P^ech^dom 

II) Phases de detection de I' instant de premiere arrivee 
d'une onde 

[0020] 

II- 1 ) Calcul de la moyenne M et de recart-type a sur 
le signal acquis en pre-declenchement sur une fe- 
netre W de dur6e w.np, ou w est une duree eiemen- 
taire modifiable a volonte. 

II-2) Formation d'un signal E c code en binaire (Fig. 
3a) obtenu en codant a 1 tout echantillon de signal 
dont le niveau en valeur absolue est superieur ou 
egal a une valeur-seuil TH1 egale a : M + n v o, et a 
0 dans le cas contraire (n-, etant un coefficient mo- 
difiable a volonte, egal a 6 dans le cas d6crit). 

II-3) Integration (Fig.3b, 4a a 4c) par sommation sur 



une fenetre glissante de largeur w.np des niveaux 
non nuls, avec affectation de la valeur obtenue dans 
chaque fenetre a rechantillon au milieu de la fenetre 
consider6e, tous les autres echantillons de cette fe- 
netre etant affectes de la valeur 0. On designe par 
S le signal forme par tous les echantillons centraux 
des fenetres glissantes successives. 

Pour diminuer les temps de calcul, ii est possi- 
ble la plupart du temps, d'augmenter le pas de glis- 
sement de la fenetre, sans pour autant alferer sen- 
siblement le r6sultat. 

II— 4) S'etant au prealable fixe pour chaque compo- 
sante un seuil relatif TH2 correspondant au rapport 
de la valeur maximale S M du signal S (Fig.4a a 4c) 
a une valeur fix6e au prealable C s1 modifiable si be- 
soin est, on detecte pour chaque composante un 
evenement representatif en identifiant le premier 
point F p depassant ce seuil TH2. Dans I'intervalle 
allant du debut du signal code somme S au point 
F P , on recherche pour chaque composante le der- 
nier point de valeur nulle LP 0 . 

II-5) On determine alors I'instant ti de premiere ar- 
rivee du signal en recherchant dans I'intervalle al- 
lant de LP 0 a la fin du signal code somme S, te pre- 
mier maximum am qui depasse le niveau TH3 = M 
+ n 2 .a, ou n 2 est comme n 1 , un entier modifiable 
(egal a 7 dans I'exemple traife). L'instant ti s'obtient 
alors en remontant le signal a partir de ce point d'un 
quart de pseudo-periode soit de np/4 points 
d'echantillonnage. Dans le cas de signaux a trois 
composantes si les valeurs determinees pour cha- 
cune des voies different, alors on applique une con- 
vention particuliere pour choisir la valeur de ti. En 
general, on prend par s6curit6 la valeur la plus bas- 
se pour ne pas manquer le premier impetus dans 
le signal. 

11-6) On peut completer les etapes pr6cedentes par 
une identification de la polarite du premier mouve- 
ment. 

[0021] En detectant avec precision I'instant de pre- 
45 miere arrivee, ti, on determine par la-meme avec une 
bonne fiabilite, la polarite du signal (positive ou negati- 
ve) aussitot apres sur chacune des voies. Si on connait 
cette polarite pour un nombre suffisant de capteurs 
ayant enregistre revenement, on peut identifier un nrfe- 
so canisme au foyer sous certaines hypotheses, notam- 
ment Identification du mecanisme bien connu dit de 
« double-couple » couramment utilise par les sp6cialis- 
tes. 

55 ill) Phases de detection de la fin et de la dur6e du signal 

[0022] L'algorithme de detection d6crit au paragraphe 
precedent peut §tre applique directement pour la detec- 
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tion de I'instant tf de fin du signal en prenant pour refe- 
rence cette fois-ci la fin du signal. 
[0023] On peut calculer la duree du signal At = tf - ti. 
On extrait le signal utile du fichier du signal initial (Fig. 
1a-1c ou 2a-2c) dans la plage allant de ti a tf avec la 5 
possibility de garder une zone de pre-declenchement 
ainsi qu'une zone de post-declenchement dont les lar- 
geurs peuvent etre proportionnelles a la pseudo-perio- 
de du signal. 

IV) Identification d'arrivee multiples distinctes au sein 
d'un meme fichier signal 

[0024] 

1 ) En pointant les differentes plages de valeurs nul- 
les du signal code somme ont peut aisement pointer 
et de ce fait comptabiliser les arrivees distinctes au 
sein d'un meme fichier signal. 

2) Une analyse en polarisation peut etre effectuee 
sur chacune des differentes arrivees pour determi- 
ner pour chacune les vecteurs propres definissant 
le triedre propre de chacune des ondes, en se re- 
ferant par exemple a celle d6crite dans le brevet FR 
2.614.997 (US 4,881 ,207) du demandeur. En com- 
parant les orientations de ces vecteurs propres, il 
est possible de detecter d'eventuelles successions 
d'ondes de compression et de cisaillement (cas ou 
les premiers vecteurs propres des deux arrivees 
consid6rees sont orthogonaux entre eux) ou bien 
d'etudier s'il s'agit de signaux correspondant a des 
multiplets (signaux presentant des signatures ana- 
logues en un meme point de reception mais acquis 
a des instants diff events). 

Avantages de la methode 

[0025] La methode selon Pinvention consiste done 
pour I'essentiel, on I'a vu, a evaluer la frequence domi- 
nante de chaque signal acquis et a utiliser la pseudo- 
p6riode correspondant a cette frequence dominante 
pour investiguer un signal integral S code en binaire et 
somme etabli a partir du signal initial brut. Par un jeu de 
quelques parametrestels que w, n 1 , n 2 etthl , ajustables 
en fonction des conditions locales d'acquisition des si- 
gnaux, on peut affiner la detection pour un jeu de don- 
nees particulier (representatif d'un site de reception de 
signaux par exemple). La methode est auto-adaptative 
car elle tient compte des signaux acquis par le biais des 
valeurs M, o, S M , F dom , np (etapes l-5, 11-1 et II-4 no- 
tamment). 

[0026] Elle permet la detection automatique de I'arri- 
vee d'une onde sismique (ou de tout autre nature) avec 
une tres bonne precision, cette precision pouvant meme 
etre amelioree automatiquement en effectuant une ana- 
lyse en polarisation a la suite de cette detection, ce qui 
est particulierement utile lorsque le debut du signal est 



noy6 dans le bruit. 

[0027] Elle contribue a la localisation automatique 
des evenements microsismiques presentant une suc- 
cession d'ondes de compression et de cisaillement. 
[0028] De meme, elle contribue a la localisation auto- 
matique des evenements sismiques ou microsismiques 
acquis sur un nombre suffisant de capteurs suffisam- 
ment espaces (4 par exemple) en permettant une trian- 
gulation automatique a partir des temps d'arrivee ainsi 
determines, le modele de vitesse etant connu par 
ailleurs. 

[0029] Elle contribue a Identification automatique 
des mecanismes au foyer lorsqu'elle est appliquee en 
un nombre suffisant de recepteurs sismiques ayant en- 
registre un evenement donne. 

[0030] La methode a ete decrite principalement dans 
le cadre duplications a I'analyse de signaux sismiques 
ou micro-sismiques. On ne sortirait pas du cadre de I'in- 
vention en I'utilisant aussi notamment en analyse de la 
parole pour le reperage de signaux successifs distincts. 



Revendications 

1. Methode d'analyse de signaux acquis provenant 
d'au moins un capteur d'ondes pour pointer auto- 
matiquement avec une grande precision sur chacun 
de ces signaux, au moins un instant significatif (ti, 
tf), caracterisee en ce qu'elle comporte : 

la determination de la pseudo-p6riode corres- 
pondant a I'inverse de la frequence dominante 
du spectre de frequence d'au moins une partie 
de chaque signal acquis, 

le codage de chaque signal acquis pour former 
un signal code en binaire par reference a une 
premiere valeur-seuil (TH1) dependant du si- 
gnal acquis, 

la formation pour chacun d'eux d'un signal in- 
tegral (S) par sommation du signal code en bi- 
naire sur une fenetre glissante, dont la largeur 
depend de la pseudo-periode, de maniere a de- 
terminer une deuxieme valeur-seuil (TH2), et le 
positionnement de I'instant significatif (ti, tf) sur 
le signal acquis par reference a un instant ou 
le signal integral correspondant atteint cette 
deuxieme valeur-seuil (TH2). 

2. Methode selon la revendication 1 , caracterisee en 
ce que le positionnement de I'instant significatif (ti, 
tf) de chaque signal acquis comporte une determi : 
nation d'un premier point de reference ou I'amplitu- 
de du signal integral correspondant (S) est egal a 
cette deuxieme valeur-seuil (TH2), la detection d'un 
premier instant intermediate (LPO) entre ce pre- 
mier point de reference et la borne la plus proche 
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du dit signal integral, la detection d'un deuxi&me 
instant intermediaire entre ce premier instant inter- 
mediate (LP 0 ) et la borne opposee du signal inte- 
gral ou un premier maximum (S M ) de ce signal in- 
tegral depasse une troisieme valeur-seuil (TH3) de- 
pendant du signal acquis, et le reperage de I'instant 
significatif recherche sur le signal acquis a un quart 
de pseudo-periode du deuxieme instant interme- 
diaire. 

3. Methode selon la revendication 1 ou 2, caracteri- 
see en ce que I'instant significatif est I'instant (ti) 
de premiere arrived de chaque signal acquis. 

4. Methode selon la revendication 2, caracterisee en 
ce que I'instant significatif est I'instant (tf) de fin de 
chaque signal acquis. 

5. Methode selon la revendication 1 , caracterisee en 
ce que Ton determine I'instant (ti) de la premiere 
arrivee du signal acquis et I'instant (tf) de fin de cha- 
que signal acquis et Ton isole sur les enregistre- 
ments des signaux regus par chaque capteur d'on- 
des, des parties d'enregistrement entre ces deux 
instants. 

6. Methode selon la revendication 5, caracterisee en 
ce que I'on isole toutes les parties d'enregistrement 
deiimitees par les dits instants significatifs, de ma- 
nure a detecter toutes les arrivees distinctes sur un 
meme capteur d'ondes. 

7. Methode selon la revendication 3, dans laquelle on 
analyse les signaux regus respectivement par une 
plurality de capteurs d'ondes relatifs a un meYne 
ev6nement g£n£rateur de signaux eiastiques, ca- 
racterisee en ce qu'elle comporte la determination 
de la polarisation de chaque signal acquis imme- 
diatement apres son instant de premiere arrivee (ti) 
et la detection des arrivees successives d'ondes de 
compression et d'ondes de cisaillement par compa- 
rison des orientations des vecteurs prop res lies 
aux polarisations respectives de ces signaux. 

8. Methode selon la revendication 3, dans laquelle on 
analyse les signaux regus respectivement par une 
plurality de capteurs d'ondes relatifs a un meme 
evenement generateur de ces signaux, caracteri- 
see en ce qu'elle comporte la determination de la 
polarisation de chaque signal acquis immediate- 
ment apres son instant de premiere arrivee (ti) et la 
detection de multiplets par comparaison des orien- 
tations des vecteurs propres N6s aux polarisations 
respectives de ces signaux. 

9. Methode selon Tune des revendications preceden- 
tes, caracterisee en ce que I'on repere des ins- 
tants significatifs distincts sur des signaux de paro- 
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1 . Verfahren zur Analyse erfasster, aus einem Wellen- 
sensor stammender Signale, urn automatisch mit 
groGer Genauigkeit an jedem dieser Signale wenig- 
stens einen signifikanten Zeitpunkt (ti, tf) zu markie- 

10 ren, dadurch gekennzeichnet, dass es umfasst: 

a. die Bestimmung der Pseudoperiode, die 
dem Kehrwert der dominanten Frequenz des 
Frequenzspektrums wenigstens eines Teils je- 

15 des erfassten Signals entspricht, 

b. die Codierung jedes erfassten Signals, urn 
unter Bezugnahme auf einen ersten, vom er- 
fassten Signal abhangigen Schwellenwert 

20 (TH1 ) ein binar codiertes Signal zu bilden, 

c. fur ein jedes hiervon die Bitdung eines Inte- 
gralsignals (8) durch Summieren des binar co- 
dierten Signals uber ein Wanderfenster, des- 

25 sen Breite von der Pseudoperiode abhangt, 

derart, dass ein zweiter Schwellenwert (TH2) 
bestimmt wird, und die Positionierung des si- 
gnifikanten Zeitpunktes (ti, tf) am erfassten Si- 
gnal unter Bezugnahme auf einen Zeitpunkt, zu 
30 dem das entsprechende Integralsignal diesen 

zweiten Schwellenwert (TH2) erreicht. 

2. Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet, dass die Positionierung des signifikanten 

35 Zeitpunkts (ti, tf) jedes erfassten Signals eine Be- 
stimmung eines ersten Bezugspunktes, bei dem die 
Amplitude des entsprechenden Integralsignals (S) 
gleich diesem zweiten Schwellenwert (TH2) ist, die 
Detektion eines ersten Zwischenzeitpunktes (LP 0 ) 
40 zwischen diesem ersten Bezugspunkt und dem 
nachsten Grenzwert dieses Integralsignals, die De- 
tektion eines zweiten Zwischenzeitpunktes zwi- 
schen diesem ersten Zwischenzeitpunkt (LP 0 ) und 
dem entgegengesetzten Grenzwert des Integralsi- 
45 gnals, wo ein erstes Maximum (S M ) dieses Integral- 
signals einen dritten, vom erfassten Signal abhan- 
gigen Schwellenwert (TH3) uberschreitet, und die 
Markierung des gesuchten signifikanten Zeitpunk- 
tes am erfassten Signal bei einem Viertel der Pseu- 
50 doperiode des zweiten Zwischenzeitpunktes um- 
fasst. 

3. Verfahren nach Anspruch 1 oder 2, dadurch ge- 
kennzeichnet, dass der signifikante Zeitpunkt der 

55 Zeitpunkt (ti) der ersten Ankunft jedes erfassten Si- 
gnals ist. 

4. Verfahren nach Anspruch 2, dadurch gekenn- 
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zeichnet, dass der signifikante Zeitpunkt der Zeit- 
punkt (tf) des Endes jedes erfassten Signals ist. 

5. Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet, dass der Zeitpunkt (ti) der ersten Ankunft 
des erfassten Signals und der Zeitpunkt (tf) des En- 
des jedes erfassten Signals bestimmt werden und 
an den Aufzeichnungen der von jedem Wellensen- 
sor empfangenen Signale Aufzeichnungsteile zwi- 
schen diesen beiden Zeitpunkten isoliert werden. 

6. Verfahren nach Anspruch 5, dadurch gekenn- 
zeichnet, dass samtliche Aufzeichnungsteile, die 
durch diese signif ikanten Zeitpunkte begrenzt sind, 
derart isoliert werden, dass samtliche unterschied- 
lichen Ankunfte an ein und demselben Wellensen- 
sor detektiert werden. 

7. Verfahren nach Anspruch 3, bei dem die Signale 
analysiert werden, die jeweils durch eine Vielzahl 
von Wellensensoren empfangen werden, die ein 
und dasselbe, elastische Signale erzeugende Er- 
eignis betreffen, dadurch gekennzeichnet, dass 
es die Bestimmung der Polung jedes erfassten Si- 
gnals unmittelbar nach seinem Zeitpunkt der ersten 
Ankunft (ti) und die Detektion der aufeinanderfol- 
genden Ankunfte von Kompressionswellen und 
Scherungswellen durch Vergleich der Orientierun- 
gen der Eigenvektoren umfasst, die mit den jewei- 
ligen Polungen dieser Signale verknupft sind. 

8. Verfahren nach Anspruch 3, bei dem die Signale 
analysiert werden, die jeweils durch eine Vielzahl 
von Wellensensoren empfangen werden, die ein 
und dasselbe, diese Signale erzeugende Ereignis 
betreffen, dadurch gekennzeichnet, dass es die 
Bestimmung der Polung jedes erfassten Signals 
unmittelbar nach seinem Zeitpunkt der ersten An- 
kunft (ti) und die Detektion von Multipletts durch 
Vergleich der Orientierungen der Eigenvektoren 
umfasst, die mit den jeweiligen Polungen dieser Si- 
gnale verknupft sind. 

9. Verfahren nach einem der vorhergehenden Anspru- 
che, dadurch gekennzeichnet, dass unterschied- 
liche signifikante Zeitpunkte an Sprachsignalen 
markiert werden. 



Claims 

1 . Method for analysis of acquired signals from at least 
one wave sensor to automatically locate with high 
precision in each of these signals at least one sig- 
nificant instant (ti, tf), characterised by the fact that 
it includes: 

determination of the pseudo-period corre- 



sponding to the inverse of the dominant fre- 
quency of the frequency spectrum of at least a 
part of each acquired signal, 

5 - coding of each acquired signal to form a binary 
coded signal by reference to a first threshold 
value (TH1) dependent on the acquired signal, 

formation for each of them of an integral signal 
10 (S) by stacking the binary coded signal in a slid- 

ing window, the width of which depends on the 
pseudo-period, so as to determine a second 
threshold value (TH2), and positioning the sig- 
nificant instant (ti, tf) of the acquired signal by 
15 reference to an instant at which the correspond- 

ing integral signal attains this second threshold 
value (TH2). 

2. Method as described in claim 1 , characterised by 
20 the fact that the positioning of the significant instant 

(ti, tf) of each acquired signal includes determina- 
tion of a first reference point at which the amplitude 
of the corresponding integral signal (S) is equal to 
this second threshold value (TH2), detection of a 

25 first intermediate instant (LPO) between this first 
reference point and the closest boundary of the said 
integral signal, detection of a second intermediate 
instant between this first intermediate instant (LP 0 ) 
and the opposite boundary of the integral signal at 

30 which a first maximum (S M ) of this integral signal 
exceeds a third threshold value (TH3) dependent 
on the acquired signal, and setting the significant 
instant sought in the acquired signal at a quarter of 
a pseudo-period from the second intermediate in- 

35 stant. 

3. Method as described in claim 1 or 2, characterised 
by the fact that the significant instant is the instant 
(ti) of first arrival of each acquired signal. 

40 

4. Method as described in claim 2, characterised by 
the fact that the significant instant is the instant (tf) 
of ending of each acquired signal. 

45 5. Method as described in claim 1 , characterised by 
the fact that the instant (ti) of first arrival of the ac- 
quired signal and the instant (tf) of ending of each 
acquired signal are determined and, in the record- 
ings of the signals received by each wave sensor, 

50 parts of recording between these two instants are 
isolated. 

6. Method as described in claim 5, characterised by 
the fact that all the parts of recording delimited by 

55 the said significant instants are isolated so as to de- 
tect all the distinct arrivals at a same wave sensor. 

7. Method as described in claim 3, in which the signals 
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respectively received by a plurality of sensors of 
waves relating to a same event generating elastic 
signals are analysed, characterised by the fact 
that it includes determination of the polarisation of 
each acquired signal immediately after its instant of 5 
first arrival (ti) and detection of the successive ar- 
rivals of compressional waves and shear waves by 
comparison of the orientations of the eigenvectors 
linked to the respective polarisations of these sig- 
nals. ™ 

8. Method as described in claim 3, in which the signals 
respectively received by a plurality of sensors of 
waves relating to a same event generating these 
signals are analysed, characterised by the fact '5 
that it includes determination of the polarisation of 
each acquired signal immediately after its instant of 
first arrival (ti) and the detection of multiplets by 
comparison of the orientations of the eigenvectors 
linked to the respective polarisations of these sig- 20 
nals. 

9. Method as described in one of the preceding claims, 
characterised by the fact that distinct significant in- 
stants are located in word signals. 25 
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